o 
o 

oo 



c/3 



C^ 



o 
o 



> 

oo 

o 
o 

-(— > 



I 

o 
o 



X 



Surface polaritons on metallic and semiconducting cylinders: 
A complex angular momentum analysis 

Stcphanc Ancey,^'|j Yves Decanini,^'[j Antoinc Folacci,^'lf| and Paul Gabrielli^'l^ 

' UMR CNRS 6134 SPE, Eqmpe Ondes et AcousUque, 

Umversite de Corse, Faculte des Sciences, Boite Postale 52, 20250 Corte, France 

^ UMR CNRS 6134 SPE, Equipe Physique Semi-Classique (et) de la Matiere Condensee, 

Umversite de Corse, Faculte des Sciences, Boite Postale 52, 20250 Corte, France 

(Dated: February 2, 2008) 

We revisit scattering of electromagnetic waves from metallic and semiconducting cylinders in the 
framework of complex angular momentum techniques. We prove that "resonant surface polariton 
modes" are generated by a unique surface wave, i.e. a surface polariton, propagating close to the 
cylinder surface. This surface polariton corresponds to a particular Regge pole of the 5'-matrix of 
the cylinder. From the associated Regge trajectory we can construct semiclassically the spectrum 
of the complex frequencies of the resonant surface polariton modes which can be considered as 
Breit-Wigner-type resonances. Furthermore, by taking into account Stokes' phenomenon, we derive 
an asymptotic expression for the position in the complex angular momentum plane of the surface 
polariton Regge pole. We then describe semiclassically the surface polariton and provide analytical 
expressions for its the dispersion relation and its damping. All these features allow us to consider the 
photon-cylinder system as a kind of artificial atom where the photon plays the role of the electron. 
Finally, we briefiy discuss the implication of our results for two-dimensional photonic crystals. 

PACS numbers: 78.20.-e, 41.20.Jb, 73.20.Mf, 42.25.Fx 



I. INTRODUCTION 

In recent years, since the work of McGurn and 
Maradudin^, there has been a growing interest in the 
study of photonic crystals containing metallic or semi- 
conducting components. The frequency-dependence of 
the dielectric function e{uj) of these materials, as well as 
the presence of a frequency range where e{uj) < 0, lead 
to new features which do not exist in conventional pho- 
tonic crystals and which are mainly linked to the presence 
of surface polaritons (SP's). With this in mind, we in- 
tend to revisit the theory of SP's propagating close to 
curved metal-dielectric or semiconductor-dielectric inter- 
faces. In this paper, we will emphasize the case of the 
circular cylindrical interface. For reviews on photonic 
crystal physics see Refs. yHI^ Among the numerous ar- 
ticles dealing with SP's in photonic crystals, we refer to 
Refs. 1,5,6 7 8 9 10 11 12 1^, i.e. to a non-exhaustive list 
of works dealing more particularly with two-dimensional 
arrays of cylinders of circular cross section for which SP's 
may be excited and which are consequently relevant to 
our study. 

SP's supported by flat metal-dielectric or 
semiconductor-dielectric interfaces can be easily de- 
scribed from a theoretical point of view and their 
properties can be obtained from rather elementary 
calculations involving homogeneous and inhomogeneous 
plane waves. It should be noted that in this context, 
an SP is clearly defined as a surface wave propagating 
close to the interface with an amplitude that decays in 
an exponential fashion perpendicularly to the interface 
and into both medium. By contrast, in the presence 
of cylindrical (this is also true for spherical) interfaces, 
the corresponding theoretical analysis is a little bit 



more complicated (see, for example, Refs. ITillTslllalTTr) : 

Bessel functions must be introduced and transcendental 
equations involving them must be numerically solved, for 
example, in order to understand resonance phenomenons 
associated with SP's. As a consequence, a clear physical 
description of scattering of electromagnetic waves from 
curved metal-dielectric or semiconductor-dielectric 
interfaces does not still exist. In the scientific literature, 
this fact sometimes leads to a semantic ambiguity with 
the expression "surface polaritons" denoting the surface 
waves propagating close to the interface as well as 
the resonant electromagnetic modes they generate. In 
order to avoid such an ambiguity and to distinguish 
between the two physical phenomenons, we shall use the 
more appropriate expression "resonant surface polariton 
modes" (RSPM's) for the latter. 

Since the sixties, mainly under the impetus of Nussen- 
zveig, asymptotic (i.e., semiclassical) techniques which 
use analytic continuation of partial-wave representations 
(Mie sums) have been developed to understand scatter- 
ing of electromagnetic waves from dielectric objects. To- 
gether these techniques form the complex angular mo- 
mentum (CAM) method. In electromagnetism, it can be 
considered as a refinement of ray optics which takes into 
account "tunneling aspects" of scattering and therefore 
includes diffractive rays associated with surface waves. 
Of course, the CAM method is an asymptotic approach 
and formally it is only valid at length scales "large" com- 
pared to the wavelength of the electromagnetic field. 

The CAM method originates from the pioneering work 
of Watsoni§ dealing with the propagation and diffraction 
of radio waves around the earth. It has since been suc- 
cessfully introduced in various domains of physics. The 
success of the CAM method is mainly due to its ability 



to provide a clear description of a given scattering prob- 
lem by extracting the physical information (linked to the 
geometrical and diffractive aspects of the scattering pro- 
cess) which is hidden in partial-wave representations and 
then to semiclassically describe resonance phenomenons. 
Here the dual structure of the S-matrix associated with 
a given scattering problem plays a crucial role. Indeed, 
the S'-matrix is a function of both the frequency u and 
the angular momentum index i. It can be analytically 
extended into the complex w-plane as well as into the 
complex £-plane (CAM plane). Its poles lying in the 
fourth quadrant of the complex w-plane are the complex 
frequencies of the resonant modes. In other words, the 
behavior of the S'-matrix in the complex w-plane permits 
us to investigate resonance phenomenons. The structure 
of the S'-matrix in the complex £-plane allows us, by 
using integration contour deformations, Cauchy's The- 
orem and asymptotic analysis, to provide a semiclassical 
description of scattering in terms of rays. In that con- 
text, the poles of the S'-matrix lying in the CAM plane 
(the so-called Regge poles) are associated with diffrac- 
tion. Of course, when a connection between these two 
faces of the S'-matrix can be established, resonance as- 
pects of scattering are then semiclassically interpreted. 
For reviews of the CAM method we refer to the mono- 
graphs of Newton ^-^ , Nussenzveig^° and Grandy^i as well 
as to references therein for various applications in quan- 
tum mechanics, nuclear physics, electromagnetism, op- 
tics, acoustics and seismology. For recent applications in 
more "exotic" contexts, we refer to Refs. 22 23 24 where 
a CAM analysis of black hole scattering and black hole 
gravitational radiation is provided and to Ref. 2^ where 
the Aharonov-Bohm effect is considered. 

As far as we know, the CAM method has never been 
used to understand scattering of electromagnetic waves 
from metallic and semiconducting objects. In fact, be- 
cause of the frequency-dependence of the dielectric func- 
tion of metals and semiconductors, the extension of the 
ideas of Nussenzveig and coworkers^*^ to such a problem 
is not quite as obvious as it seems at first sight. In this 
article, we make some steps in that direction but with a 
rather modest goal. Indeed, we only consider the scatter- 
ing of TE waves {H polarization) by a metallic or semi- 
conducting circular cylinder surrounded by a dielectric 
medium. We limit our study to that case because SP's 
are not excited in the ^^-polarization configuration. From 
the S'-matrix of the cylinder, using CAM techniques, 
we develop a semiclassical description of the scattering 
aspects linked to SP's. More precisely, we prove that 
RSPM's are generated by a unique SP propagating close 
to the cylinder surface. This surface wave is associated 
with a particular Regge pole of the S'-matrix of the cylin- 
der. From the corresponding Regge trajectory, i.e. from 
the curve traced out in the CAM plane by this Regge 
pole as a function of the frequency, we can construct 
semiclassically the spectrum of the complex frequencies 
of RSPM's which can be considered as Breit-Wigner-type 
resonances. Furthermore, by carefully taking into ac- 



count Stokes phenomenon, we derive an asymptotic ex- 
pression for the position of the SP Regge pole in the CAM 
plane and then, we can describe semiclassically the SP. In 
some sense, our results allow us to consider the photon- 
cylinder system as an artificial atom: RSPM's are long- 
lived quasibound states for this atom while the trajectory 
of the SP which generates them and which is supported 
by the cylinder surface is a Bohr-Sommerfeld-type orbit. 
Our paper is organized as follows. In Section 2, we in- 
troduce our notations and we construct the S'-matrix of 
the system. We then discuss the resonant aspects of our 
problem. In Section 3, by using CAM techniques, we es- 
tablish the connection between the SP propagating close 
to the surface cylinder and the associated RSPM's. In 
Section 4, we describe semiclassically the SP by provid- 
ing analytic expressions for its dispersion relation and its 
damping. We then deduce analytic approximations for 
the excitation frequencies of RSPM's. Finally, in Section 
5, we conclude our article by considering some possible 
extensions of our work and by briefly discussing the im- 
plication of our results in the context of two-dimensional 
photonic crystal physics. 



II. EXACT S'-MATRIX AND SCATTERING 
RESONANCES 

From now on, we consider the scattering of an electro- 
magnetic wave by a metallic or semiconducting circular 
cylinder with a frequency-dependent dielectric function 
ec{uj) which is embedded in a host medium of infinite 
extent with constant dielectric function eh (region I). In 
the usual cylindrical polar coordinate system (p, 9, z) the 
cylinder occupies a region corresponding to the range 
< p < a (region II). We also assume that the mag- 
netic field H is parallel to the axis of the cylinder {H 
polarization) and we choose to treat our problem in a 
two-dimensional setting, ignoring the z coordinate. Fur- 
thermore, in the following, we implicitly assume the time 
dependence exp(— ia;i) for the magnetic field and we shall 
sometimes use the wave numbers 

k\Lu) = (^^)^, and fc"^^) = (^) v/^IR- (1) 

Here c is the velocity of light in vacuum. 

As far as the dielectric function of the cylinder is con- 
cerned, we assume it presents a Drude-like behavioriSSi^ 



ec(w) = Coo 1 1 



or the ionic crystal behavioriSSiSi 



ec{uj) = Cc 



(2) 



(3) 



In both cases, e^o is the high-frequency limit of the di- 
electric function. In Eq. |(2Jl, tjp is the plasma frequency. 



In Eq. Q , ijJT and lol respectively denote the transverse- 
optical-phonon frequency and the longitudinal-optical- 
phonon frequency. In the first case, SP's follow from the 
coupling of the electromagnetic wave with the plasma 
wave and are usually called surface plasmon polari- 
tons. In the second one, SP's follow from the coupling 
of the electromagnetic wave with the longitudinal and 
transverse acoustic waves and are usually called surface 
phonon polaritons. Eq. (0) can be used to describe the 
dielectric behavior of certain metals and semiconductors 
(Si, Ge, InSb, . . . ) while Eq. |(2Jl can be used to investi- 
gate the optical properties of other semiconductors such 
as GaAs. 

From Maxwell's equations it is easy to show that the 
z-component of the magnetic field satisfies the Helmholtz 
equation 



0%.(., 



if^^(x) =0 for < p < a. 



(4a) 



A, 






Hlixj^O forp>a, (4b) 



where x = (p, 9) and with the Laplacian Ax given, in the 
polar coordinate system, by 



A. = 



9p2 



ld_ 
pdp 






(5) 



Furthermore, from the continuity of the tangential com- 
ponents of the electric and magnetic fields (i.e., ol Eg and 
Hz) at the interface between regions I and II, it can be 
shown that the z-component of the magnetic field satis- 
fies for Q<e <2tt 



Hl(p^a,e)^Hl\p^a,0), 



1 dHl 
e/i dp 



{p^a,9) 



1 dHl 



ec{uj) dp 



(6a) 
(p^a.d). (6b) 



We are first interested in the construction of the S- 
matrix for the cylinder. Because of the cylindrical sym- 
metry of the scatterer, the S'-matrix is diagonal and its 
elements Se£' are given by Se£' — Sg Sw. For a given 
angular momentum index £ G Z, the coefficient Se is 
obtained from the partial wave (i?z)f solution of the fol- 
lowing problem (here we extend, mutatis mutan dis, the 
quantum mechanical approach developed in Ref . |28|) : 

(i) {Hz)f, satisfies the Helmholtz equation Q), 

(u) {Hz)i satisfies the boundary conditions ©, 

(iii) at large distance, {Hz)^ has the asymptotic behavior 



{Hz),{p.e) 



,-i(fcV-^7r/2-7r/4) 






Outside the cylinder (region I), the solution of Q) is ex- 
pressible in terms of Bessel functions (see Ref. 1220 as a lin- 
ear combination of Ji{k^p)e^^^ and H\ ' {k^p)e^^^ . Inside 
the cylinder (region II), it is proportional to Ji{k^^ p)e^^^ . 
As a consequence, the partial wave {Hz)i solution of (i) 
and (ii) can be obtained exactly. Then, by using the stan- 
dard asymptotic behavior of Hankel functions for a; ^ cx3 
(see Ref. Hi) 

iJ,('^(x)^V27(^e^*--^^/2-/4), (7a) 

^f'(^)^V2/(7ra;)e-'("-^'^/2-"/4), (7b) 

we find from (iii) the expression of the diagonal elements 
Sf of the S'-matrix. We have 



S({uj) 



1 



Di{lo) 



(8) 



where D\ {u}) and Df{uj) are two 2x2 determinants 
which are explicitly given by 

D^^\lu) = k^\Lu)J'^ik\u)a)Mk^^uj)a) 

-k^{uo)Ji{k^{uo)a)J'^{k^^{uj)a), (9a) 

D^iu) = k^^{uo)Hf^' {k^{uj)a)Ji{k^^{u)a) 

(9b) 



~k\uj)H)'\k\uj)a)J'^{k^^{u)a). 

The unitarity of the S-TnaXvh^, which expresses the en- 
ergy conservation, can be easily verified from (jSJ and © 
by using elementary properties of Bessel functions. The 
reciprocity property ^^, which is associated with time- 
reversal invariance, is also satisfied because S^ is an even 
function of L 

The 5-niatrix is of fundamental importance because 
it contains all the information about the scattering pro- 
cess. Its components appear in the Green functions of 
the problem, in the scattered field when a plane wave 
excites the cylinder as well as in both the scattering am- 
plitude and the total scattering cross section. As far as 
the scattering by a plane wave propagating along the x 
axis is concerned, the total magnetic field in region I is 
given by 

Hl{p,e)^e''''P^°^'^ 

+ CXD 

+ Ey*' [s,-i]Hf\k'p)cosm. (10) 

1=0 

Here 7^ is the Neumann factor (70 = 1 and for t ^ Q, 
7£ = 2). The scattering amplitude /(w, 9) is defined from 
the asymptotic behavior of the total magnetic field by 



^1(P,< 



e'^'p""""^ + f{u,9y 






(11) 



By using the asymptotic behavior (|7a|l in Eq. (|10|l . we 
can write 



f{^,0) 



1 



2iTrk^{uj) 



+ 00 

■E 



ji[Si{u;)-l]cos{£d). (12) 
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FIG. 1: a) Total cross section ar- b) Scattering resonances in 
the complex tja/c- plane. €c{io) has the Drude type behavior 
with 600 = 1 and uj^a/c = 2-n while e^ = 1. 
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FIG. 2: a) Total cross section or. b) Scattering resonances in 
the complex oja/c-plane. ec(cj) has the ionic crystal behavior 
with too = 1, loto./c — 2-K and u}lo./c = Stt while e/i = 1. 



Then, the total scattering cross section per unit length 
of the cylinder can be obtained by using the optical 
theoremiS: 



Mu^) = \/fc^I^ (e-^^/V(c^, 0^0)). (13) 

In Figs. ^ and 12^, we present two examples of total 
cross section. They are both plotted as functions of the 
reduced frequency loaj c. In Fig. ^, the cylinder is em- 
bedded in vacuum (e/j — 1) and its dielectric function is 
given by Q with try^ — \ and uj^ajc = 27r. In Fig.O^, the 
cylinder is embedded in vacuum (e^ — 1) and its dielec- 



tric function is given by Q with Coo = 1, tOTd/c = 27r and 
LUi^a/c — Stt. Even if we restrict ourselves to those par- 
ticular configurations, it should be noted that the results 
emphasized numerically and that we shall now discuss 
are very general. On the two figures, rapid variations 
of sharp characteristic shapes can be observed. This 
strongly fluctuating behavior is due to scattering reso- 
nances. These resonances are the poles of the S'-matrix 
lying in the fourth quadrant of the complex w-plane and 
they are determined by solving 



Di{uj) = for ^ G N. 



(14) 



= ^en - «r£p/2 



The solutions of (|14|l are denoted by 

where uigj > and Tip > 0, the index p permitting us to 
distinguish between the different roots of H14|l for a given 
i. In the immediate neighborhood of the resonance LOip, 
Si{u}) has the Breit-Wigner form, i.e., is proportional to 



rfp/2 



,(0) 



iTip/2 



(15) 



As a consequence, when a pole of the S'-matrix is suffi- 
ciently close to the real axis in the complex w-plane, it 
has an appreciable influence on the scattering amplitude 
and therefore on the total cross section. Of course, if 
a pole is very close to this axis, the corresponding peak 
is too sharp to be observed on the total cross section. 
In Figs.^ and 135, resonances are exhibited for the two 
configurations previously considered. A one-to-one cor- 
respondence between the peaks of ctt and the resonances 
near the real cja/c-axis can be clearly observed in certain 
frequency ranges. 

More precisely and more generally, for the dielectric 
function |2Jl there exists in the frequency range cu < ujp 
where ec{uj) < a family of resonances close to the real 
axis of the complex w-plane which converges, for large £, 
to the limiting frequency Ws satisfying 



and given by 



,(^5) +6,, = 



^/T+Thjec 



(16) 



(17) 



For the dielectric function |(3l, the same type of behav- 
ior occurs but in two different frequency ranges: in the 
region to < ujt (here ec(ct;) > 0) with an accumulation of 
resonances at the pole ut of the dielectric function and 
in the region lut < lu < lul (here ec(w) < 0) with an 
accumulation of resonances at the limiting frequency los 
still satisfying (|16|l but which is now given by 



IujI + {eh/e^)uj^ 



1 + e/i/cc 



(18) 



We must keep in mind that in the scattering of a H- 
polarized photon with frequency ujg° , a decaying state 



of the photon-cylinder system is formed. It has a finite 
hfetime proportional to 1/Tip. The resonant states whose 
complex frequencies belong to one of the families previ- 
ously described are therefore long-lived states. Because of 
these particular quasibound states, the photon-cylinder 
system behaves as a kind of artificial atom for which the 
photon plays the usual role of the electron. We shall 
come back to this point of view in the next two sections. 
From now on, we shall more particularly focus our at- 
tention on the physical interpretation of the long-lived 
resonant states whose excitation frequencies belong to 
frequency ranges in which ec(w) < 0. We shall prove 
that these states are generated by a SP propagating close 
to the cylinder surface and for this reason we call them 
RSPM's. For such states, the artificial-atom point of 
view can be pushed farther as we shall show in Sections 
3 and 4. For the long-lived resonant states whose excita- 
tion frequencies belong to the frequency range in which 
ec{uj) > (the so-called bulk polariton states), we are 
not able to provide a similar analysis. This is not very 
serious as they do not have, in photonic crystal physics, 
the importance of RSPM's. 



III. SEMICLASSICAL ANALYSIS: FROM THE 

SP REGGE POLE TO THE COMPLEX 

FREQUENCIES OF RSPM'S 

Using the CAM method, we can provide a physical pic- 
ture of the scattering process in term of diffraction by sur- 
face waves and more particularly a physical explanation 
of the excitation mechanism of RSPM's valid for "high 
frequencies" . By means of a Watson transformationiSi 
applied to the scattering amplitude 112|) . we can write 



fito,0) 



27rfci(w) jc 



(^aH-1) 
sinTrA 



cos [A(7r - 61)] dX. 

(19) 
Here C is the integration contour in the complex A- 
planei^ illustrated in Fig. O and which encircles the real 
axis in the clockwise sense. In Eq. (|19ll . V which stands 
for Cauchy's principal value at the origin is used in order 
to reproduce the Neumann factor. The Watson trans- 
formation has permitted us to replace the ordinary an- 
gular momentum £ by the complex angular momentum 
A. S\{uj) is now an analytic extension of Si{(jj) into the 
complex A-plane which is regular in the vicinity of the 
positive real A axis. Using Cauchy's theorem and by not- 
ing that inside the contour C, the only singularities of the 
integrand in (|19|) are the integers, we can easily recover 
(II21 from (PJ. 

We can then deform the path of integration in H19|) 
taking into account the possible singularities. The only 
singularities that are encountered are the poles of the 
S-matrix lying in the CAM plane. They are known as 
Regge polesiSiSS and are determined by solving 



Im A 



Re A 



FIG. 3: The Watson integration contour. 



Figs. 0] and |S1 exhibit the distribution of Regge poles for a 
cylinder embedded in vacuum. We still consider the two 
configurations previously studied. We do not display the 
Regge pole distributions for other configurations (i.e. for 
other values of the parameters Coo, e/i, wp, ujt or lul) 
because they are not really different from those of Figs.^ 
and O In fact, all these Regge pole distributions are 
rather similar to the distributions associated with the 
dielectric objects usually studied^*'. However, in the fre- 
quency range where e{uj) < 0, something new occurs: it 
exists a particular Regge pole lying in the first quadrant 
of the A-plane and very close to the real axis. It is not 
present for ordinary dielectric objects. As we shall see 
below, this new Regge pole is associated with a SP orbit- 
ing around metallic or semiconducting cylinders. From 
now on, we shall denote it by Agp(cL;). 

By Cauchy's theorem we can then extract from H19() 
the contribution of a residue series over Regge poles. In 
fact, we limit our study to the contribution of Agp(cL;) 
which is given by 



/sp('^,^) = 



27r 



rsp(cj) 



ik^{uj) sin [7rAsp(w)] 



cos [Asp(i-ij)(7r — 



Here rsp{u}) ~ residue {S\{uj))-)^^ 



A=Asp(w)' 



Dx{uj) = for uj>0. 



(20) 



(21) 
Of course, / dif- 
fers from /sp by a smooth background integral and by 
the contributions of all other Regge poles. We are not 
interested by these contributions which do not play any 
role in the excitation of RSPM's. We think that these 
contributions could be studied, mutatis mutandis, in the 
framework of CAM techniques developed by Nussenzveig 
and coworkers^. By using 
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FIG. 4: Regge poles in the complex angular momentum 
plane. ec{u}) has the Drude type behavior with eoo = 1 and 
ujpa/c — 2-K while e^ = 1. The distribution corresponds to 

uoa/c = 4. 



FIG. 5: Regge poles in the complex angular momentum 
plane. ec{i-o) has the ionic crystal behavior with Eoo = 1, 
ujTa/c — 2-K and LOj^ajc — Stt while e^ = 1. The distribution 
corresponds to coa/c — 7.8. 



which is true if Im A > 0, we can write 



/sp(c^,0)-- 




rsp(w) 



iAsp(c^)(e+2m7r) 



^iAsp(a))(27r-e+2m7r) 



(22) 



In Eq. (|22|l . exponential terms correspond to surface wave 
contributions. Because the Regge pole Asp(w) hes in the 
first quadrant of the CAM plane, exp[iAsp(a;)(6')] (resp. 
exp[iAsp('^)(27r— 6*)]) corresponds to a surface wave prop- 
agating counterclockwise (resp. clockwise) around the 
cylinder and Re Asp('^) represents its azimuthal prop- 
agation constant while Im Asp(aj) is its damping con- 
stant. The exponential decay exp[— Im Asp(w)^] (resp. 
exp[— Im Xsp{lo){2t: — 6)]) is due to continual re-radiation 
of energy. Moreover, in (|22|l . the sum over m takes into 
account the multiple circumnavigations of the surface 
waves around the cylinder as well as the associated ra- 
diation damping. The Regge pole Asp is very close to 
the real axis in the complex A-plane. It then corresponds 
to a surface wave which is slightly attenuated during its 
propagation and which contributes significantly to the 
scattering process and to the resonance mechanism. 

As cj varies, the Regge pole Asp((^) describes a Regge 
trajectory-^^ in the CAM plane. In Figs.Eland[71 we have 
displayed the Regge trajectories of SP's for the two con- 
figurations previously studied. It should be noted that 
a.s u! ^ cus, the real part of the SP Regge pole increases 
indefinitely while its imaginary part vanishes. For other 
configurations (i.e. for other values of the parameters eoo, 
eh, Wp, (jJt or ujl), we have verified that the SP Regge 
pole behavior is very similar. 

The resonant behavior of the cylinder-photon system 



can now be understood in terms of the Regge trajectory 
of the SP. When the quantity Re Asp(ti^) coincides with 
an integer, a resonance occurs. Indeed, it is produced by 
a constructive interference between the different compo- 
nents of the surface wave, each component corresponding 
to a different number of circumnavigations. Resonance 
excitation frequencies u!ggp are therefore obtained from 
the Bohr-Sommerfeld type quantization condition 



Re 



Asp U 



,{o) 
'fSP 



£ = 0,1,2, 



(23) 



(o) 



By assuming that ui is in the neighborhood of w^sp ^^^ 
using Re Asp(w) ^ Im Asp(w) (which can be numeri- 
cally verified, except for low frequencies), we can expand 

Asp(t^) in a Taylor series about i^) gp, and obtain 



Xspiiu)^e + 



dReAsp(w) 



div 



,(o) 



{u 






i ImAgp(w)g'p) -f 



(24) 



Then, by replacing H24|) in the term cos [7rAsp(ti')] of 121|l . 
we show that fsp{i-o, 9) presents a resonant behavior given 
by the Breit-Wigner formula H15|) with 



fSP 



Im Asp(w) 



d Re Asp(w)/dw 



(25) 



Eqs. (|23|) and H25|) are semiclassical formulas which per- 
mit us to determine the location of the resonances from 
the Regge trajectory of Asp. 

Tables III and ITTl present samples of complex frequencies 
of RSPM's for the two configurations previously consid- 
ered. They are calculated from the semiclassical formulas 
(I23|l and H25|) by using the Regge trajectories numerically 
determined by solving (|20|) (see Figs. [Sj and [7|). A com- 
parison between the "exact" and the semiclassical spectra 




coa/c 



FIG. 6: Regge trajectory for the SP Regge pole. ec{uj) has 
the Drude type behavior with Eoo = 1 and Upa/c = 2n while 
eii — 1. As Loa/c -^ LUga/c, the real part of the SP Regge pole 
increases indefinitely while its imaginary part vanishes. 




FIG. 7: Regge trajectory for the SP Regge pole. ec{to) has 
the ionic crystal behavior with too = 1, uira/c = 2n and 
LULd/c = Stt while e^ = 1- As uia/c — > ujaa/c, the real part of 
the SP Regge pole increases indefinitely while its imaginary 
part vanishes. 



TABLE I: The first complex frequencies of RSPM's. €c{uj) 
has the Drude type behavior with too = 1 and uOpa/c — 2-k 
while eh = 1. 





Exact 


Exact 


Semiclassical 


Semiclassical 


I 




rffip/2 


(o) 


r,sp/2 


1 


0.570278 


0.642122 


0.665828 


0.582435 


2 


1.53524 


0.633630 


1.49245 


0.602114 


3 


2.35038 


0.475286 


2.25693 


0.476399 


4 


2.97766 


0.274602 


2.90439 


0.289791 


5 


3.42645 


0.119226 


3.39638 


0.127238 


6 


3.73632 


0.036765 


3.73016 


0.037645 


7 


3.94064 


0.007494 


3.94009 


0.007435 


8 


4.07049 


0.000999 


4.07047 


0.000981 


9 


4.15483 


0.000093 


4.15476 


0.000094 


10 


4.21272 


0.000006 


4.21266 


0.000006 



TABLE II: The first complex frequencies of RSPM's. f.c(uj) 
has the ionic crystal behavior with Ecx) = 1, uira/c = 27r and 
lolcl/c = 37r while e^ = 1. 



Exact 



Exact 

r^sp/2 



Semiclassical 

,{o) 



LOi 



Semiclassical 

r«sp/2 



7 


6.459387 


0.4345112 


6.283821 


0.0837525 


8 


6.809183 


0.2419456 


6.698432 


0.3004512 


9 


7.091176 


0.1055079 


7.056186 


0.1188560 


10 


7.311233 


0.0353546 


7.304271 


0.0371596 


11 


7.470835 


0.0088275 


7.470046 


0.0088537 


12 


7.581710 


0.0016441 


7.581652 


0.0016417 


13 


7.659584 


0.0002366 


7.659579 


0.0002392 


14 


7.716557 


0.0000273 


7.716543 


0.0000274 


15 


7.759919 


0.0000026 


7.759919 


0.0000027 



IV. SEMICLASSICAL ANALYSIS: 

ASYMPTOTICS FOR THE SP AND PHYSICAL 

DESCRIPTION 

A deeper understanding of the SP behavior can be ob- 
tained by solving perturbatively Eq. l(^ for A = Asp- 
We first replace the Bessel function J\{z) by the modi- 
fied Bessel function Ix{z) (see Ref. [297 in order to take 
into account the fact that ec(w) < 0. Eq. if^ reduces to 



shows a very good agreement, except for "low" frequen- 
cies. We have also performed the corresponding calcula- 
tions for other configurations with Ch ^ ^ and Coo 7^ 1- 
The agreement seems even better. Furthermore, inserted 
into the semiclassical formulas (|23|l and H25|) . the behav- 
ior of Regge trajectories near the limiting frequencies uJs 
easily explains the existence of the families of resonances 
close to the real axis of the complex w-plane which con- 
verges for large i to the limiting frequency ojs- In con- 
clusion, we have established a connection between the 
complex frequencies of RSPM's and a particular surface 
wave, the so-called SP, described by a particular Regge 
pole of the S'-matrix and which orbits around the cylin- 
der. 



1 HiZiV^a^/c) _ 1 r^^^i^-e,{Lo)aLu/c) 

(26) 
In the right-hand side of H2t)|l . we can use the uni- 
form asymptotic expansions oi I\{z) for large orders (see 
Ref.l^, i.e. 



h{zy 



1 



1 



where 

Fxiz) 



^(A2+z2)l/4 

= (A2 + z2)i/2 + Aln(' 



:.Fxiz)/2 



A+(A2+z2)l/2 



(27) 



(28) 



By assuming |Asp| ^ ^J--tJ^Lo)aui j c^ we then obtain 

1 ^Asp(v/-^^a'^/c) 



^-ec(t^) /Asp(V-ec(w)aw/c) 



[A 



•SP 



■ec{w)(aw/c)2]^ ^ 



c(a;)(aa;/c) 



(29) 



In the left-hand side of (|26|l . the relative positions of Asp 
and ^JThau I c in the A-complex plane (see Fig. llO|) permit 

us to employ the Debye asymptotic expansion of _ff| (z) 
in the form22i^ 



fff^(z)--iA(A,z)e-"(^'^) 



where 

A(A,.)=(^) {\--z^)-^l\ 

'z + (A2 



(30) 



(31a) 



i(A,z) = (A2-z2)V2_Aln(': 



^2)1/2 



By assuming |Asp| 3> ^Jchauj/c, we then deduce 

1 H^^li^uaoj/c) [\l^~eh{aoj/cf]"^ 



^<Uv^«^/c) 



e,,(aw/c) 



Eq. I|26|l can now be solved and we easily find 



ASP(^) - (-) \^^y 



(31b) 



(32) 



(33) 



We have obtained an asymptotic expansion for Agp or 
more exactly for the real part of Asp- Indeed, the right- 
hand side of (|33|l is purely real. The perturbative resolu- 
tion of Eq. (|26|l did not permit us to extracted the small 
imaginary part of Asp- Of course, it would be possible 
to improve (|33|) by taking into account higher orders in 
the asymptotic expansions (|27|l and (|3L)|I . But that does 
not seem necessary. First, this does not provide an imag- 
inary part for Asp. In fact, as we shall see below, this 
term corresponds to an exponentially small contribution 
which lies beyond all orders in perturbation theory. Fur- 
thermore, we have numerically tested the formula (|33|l for 
various values of the parameters Eoo, £fn Wp, lot or lul. 
In all cases, it provides a rather good approximation for 
Re Asp (see, for example. Figs. |S1 and 1^1 for the two con- 
figurations previously studied). It should be noted that 
it also predicts the divergence of Re Asp for to ^ uJs- 

We can now insert (|33|) into (|22|l . The contributions 
corresponding to the SP propagating counterclockwise 
and clockwise are then given by 

exp [i{±\sp{uj)9 ~ cut)] = exp [i{±ksp{uj)a0 - ut)] (34) 

with 



fcsp(-) = Q ^ 



eh + ec{uj)' 



(35) 




6 (0 a/c 



FIG. 8: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. ec{io) has the Drude 
type behavior with e^o ~ 1 and uj^ajc — 2tt while e^ = 1. 




FIG. 9: Regge trajectory for the SP Regge pole. Comparison 
between exact and asymptotic theories. ec(a)) has the ionic 
crystal behavior with too = 1, uJTa/c = 2tt and lulcl/c = 3n 
while e^ = 1. 



Here we have reintroduced the time dependence 
exp(— iwt) in order to clarify the physical interpretation. 
From (|34|) and by noting that a dO represents the length 
element on the cylinder surface, it now appears that 
the SP propagation is supported by the cylinder surface 
which thus plays the role of a Bohr-Sommerfeld-type or- 
bit and that (|35|l can be considered as the SP dispersion 
relation. This relation could permit us to derive analyt- 
ically the phase velocity Vp = uj/kQp[ijj) as well as the 
group velocity Vg = duj/dksp{uj) of the SP. 

It should be also noted that the dispersion relation l|35|) 
is in fact the usual dispersion relation of a SP supported 
by a flat metal-dielectric or semiconductor-dielectric in- 
terface (see, for example Ref. l3^ . We have recovered 
the same dispersion relation because we have limited the 
perturbative resolution of Eq. H2t)|l to the lowest order. 
By taking into account higher orders in the asymptotic 



expansions (|27|l and (|30() , we could obtain corrections for 
(|33|l and l|35|) vanishing for a — > c», i.e. in the hmit of 
large radius. 

By inserting the expression H33I) for Asp(w) into the 
Bohr-Sommerfeld quantization condition H23|l . we can de- 



rive approximations for the resonance excitation frequen- 
cies w)gp. If the dielectric function of the cylinder is given 
by ©I we obtain for the reduced frequencies 



,(o) 



1 

71 



(h 



f-h(-c 



e- 



m 



e/i 



e^eo 



1 2 



^2 



e/i 



£2 



1/2 



This analytic formula provides accurate results for "large" values of £. For 



(36) 



= 3, the error is around 13 % and it 

X). For 



becomes less than 1 % for £ > 7. Furthermore, it predicts the convergence to uJs = ^p/-\/l + Qi/^oo when £ 
the dielectric function (|5Jl, we obtain 



U} 



esp 



1 

V2 



V c 



c 



e/j 



e/iCo 



e/i 



e/jeo 



^2 



e/jec 



(^)'+«(T)^ 



^2 



1/2 



(37) 



r 



ii Im A 



Reo- >0 



Re a <0 



SP 
O 



Re A 



Stokes line 
Im a =0 

Principal curve 
Rece =0 



Re a <0 



FIG. 10: The relative positions, in the A-complex plane, of 
the Regge pole Asp and the reduced frequency /? = y/ehua/c. 



This analytic formula provides accurate results for 
"large" values of i. For i — 7, the error is around 
3 % and it becomes less than 1 % for £ > 10. Fur- 
thermore, it predicts the accumulation of resonances at 
UJs = ^JiiJi -I- (e,j/eoo)c^r)/(l + e/i/eoo) for ^ -> oo. 

To conclude this section, we shall come back to the 
damping constant Im Asp(w) of the SP. Numerically (see 
Figs. El and O we have shown that this term is small 
but we cannot consider that it vanishes as we previously 
found. In fact, this term corresponds to an exponentially 
small contribution which lies beyond all orders of the 
asymptotic expansion (|30|l and which can be captured by 



carefully taking into account Stokes' phenomenon^i^. 
(For modern aspects of asymptotics beyond all orders and 
of Stokes' phenomenon, we refer to a beautiful article by 
ggj.j.y35 g^g .^gjj g^g -(-Q Ref s . .3 6^3 7.3 S^ ■ Instead of (I30|l . we 

shall use the Debye asymptotic expansion of H\ {z) in 
the form 



h[^\z)^ - iA{X, z)e-"(^^") (1 + . . . ) 

+S [a(A, z)] A{\, z){l + ... )e"(^^^). (38) 



In the right-hand side of 138|l . the first term is the usual 
Debye asymptotic expansion truncated near its least 
term. The second one is obtained by decoding the diver- 
gent tail of that asymptotic expansion. This can be done 
(see Refs. 35 38) by Borel summation after exploiting a 
resurgence formula discovered by Dingle'^^. In the region 
of the A-complex plane where the Regge pole Asp lies 
(see Fig. llU|l . we have Re a < 0. As a consequence, the 
first term of the right-hand side of 138|) is the dominant 
contribution while the second one is a subdominant term 
which can be forgotten when |A| — > oo. That is what we 
did previously by using (|30|l . The Stokes multiplier func- 
tion S [a(A, z)] is a complicated function involving the 
exponential integral function Ei. It goes continuously 
from to 1 at the crossing of the Stokes line Im a = 
emerging from the turning point z = (3 (see Fig. llOf) . Be- 
low the Stokes line, it rapidly vanishes. On the Stokes 
line, it is equal to 1/2 and above the Stokes line it rapidly 
becomes equal to 1. It thus describes the rapid but con- 
tinuous birth of the subdominant contribution near the 
Stokes line. 
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From H38() we can now write 

1 H'^^li^nauj/c) [\l^ - en{aio/cf]^'^ 



r(i) 



hH):>^{^aLo/c) 



ehiau/c) 



(l - 2iS [a(Asp, A/eft^aw/c)] e 



2a(Asp ,y/ehaLj/c) 



(39) 



instead of (EH). By using (EU and (EH, Eq. (ESJ can be 
solved approximately and we find 



Re 



Asp(c.) - (^) J^f^ 



(40a) 



ImAsp(w) - 2 (— ) P(cj)5(a;)e2"(") (40b) 



where 



Piu) 



eUii^) 



(4 - e2(w)) y^eh€c{uj){€h+ec{uj)) 



, (41a) 



a(c^) ^ a l{u;a/c)J "'''^''.\ , ^auj / c] , (41b) 



Eh + ec{u)) ' 



S{uj) « S 



a {uja/c)\ — ^/ , N : V^au/c 



(41c) 



We can see easily that the imaginary part (j40bp of Asp 
vanishes for u; = cOs bs well as in the large radius limit 
a — > oo, i.e. in the flat interface limit. Furthermore, we 
have numerically studied Ij40b|) for "high" values of w, i.e. 
when the Regge pole Asp is very close to the Stokes line. 
In that case, by giving to the Stokes multiplier function 
5 the value 1/2, we have checked that (j40b|) provides 
accurate results for the imaginary part of Asp . We there- 
fore consider we have succeeded in providing an analytic 
formula for the Regge pole of the SP. However, the ex- 
pression (j40b|) is a rather complicated function of w. As 
a consequence, its use in the semiclassical formula (P5|) 
is unfortunately not very interesting. In short, we think 
that (|40bp is especially interesting for the qualitative de- 
scription of the SP damping it provides. 



V. CONCLUSION AND PERSPECTIVES 

In this article we have introduced the CAM method in 
the context of scattering of electromagnetic waves from 
metallic and semiconducting cylinders. This allows us to 
provide a physical explanation for the excitation mech- 
anism of RSPM's as well as a simple mathematical de- 
scription of the surface wave (i.e. the SP) that generates 
them. It should be noted that our results are not limited 
to metals and semiconductors. Under simple assump- 
tions, they are also valid, mutatis mutandis, for more 
general materials. Indeed, in a frequency range where 



the dielectric function of a material presents a dominant 
simple pole lUq, it is always possible to write^ 



ec{i^) 



2ujoR 



iJl — Lj(ljJ + 27) 



(42) 



Here Wo is the resonance frequency of the material in 
the frequency region considered, 7 denotes the associated 
damping term, too stands for the high-frequency limit of 
the dielectric function and we assume that the coefficient 
R is positive. In the absence of dissipation (7 w 0) and if 
the zero of tJyuS) which is given by y/oj^+%j^Jl/€^ lies 
in the validity range of Eq. (|42|l , there exists a SP which 
can be described by the theory developed in Sections 3 
and 4. 

In parallel with the semiclassical analysis of SP's on 
metallic or semiconducting cylinders, we have developed 
a new picture of the photon-cylinder system: it can be 
viewed as an artificial atom for which the photon plays 
the role of an electron. RSPM's are long-lived qua- 
sibound states for this atom, the associated complex 
frequencies are Breit-Wigner-type resonances while the 
trajectory of the SP which generates them is a Bohr- 
Sommerfeld-type orbit. Furthermore, the imaginary part 
of a given RSPM complex frequency corresponds to an 
exponentially small term which lies beyond all orders in 
perturbation theory. As a consequence, as their exci- 
tation frequency increases, RSPM's gain very long life- 
times, i.e. they behave like bound states. 

With in mind applications in photonic crystal physics, 
our work could be naturally extended in various direc- 
tions including i) scattering by cylinders with metallic 
or semiconducting coating, or more generally, multilay- 
ered structures, ii) scattering by metallic or semiconduct- 
ing spheres and last but not least iii) scattering by ob- 
jects fabricated from left-handed materials. It would be 
also interesting to provide a complete (i.e., not limited 
to SP's) semiclassical description of scattering of elec- 
tromagnetic waves from metallic and semiconducting ob- 
jects in the framework of CAM techniques by extending 
the ideas of Nussenzveig and coworkersSS. But at first 
sight that seems to be a formidable task. 

In recent papers dealing with photonic band struc- 
ture of two-dimensional photonic crystals fabricated from 
metallic or semiconducting cylinders arrayed in a square 
lattice, a striking feature has been notedii^iSi^iiSiiiiiSii^, 
namely the existence of flat bands (i.e. dispersionless 
bands) in the frequency range in which e{uj) < 0. This 
result, which only exists for H polarization, is of course 
linked to the excitation of RSPA4's. More precisely, it is 
due to the localization of the photon which is trapped 
on the Bohr-Sommerfeld orbit. Of course, this analy- 
sis is rather over-simplified. In fact, it is necessary to 
understand up to what point single-cylinder resonant as- 
pects are related to "resonant" aspects of the full pho- 
tonic crystal. Recently, Ito and Sakodc^iii have consid- 
ered this problem by developing a physically intuitive but 
appealing analysis: they regard the RSPM's of an iso- 
lated cylinder as atomic orbitals and they describe their 
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effects into the photonic crystal in the context of the "hn- 
ear combination of atomic orbitals (LCAO) theory" . Ito 
and Sakoda do not use the terminology "artificial atom" 
to describe the photon-cylinder system but this picture is 
implicitly present in their work and the point of view we 
develop in the present article strengthens their analysis. 
Of course, it should be interesting to provide a more 
rigorous interpretation of the existence of the flat bands. 
With this aim in view, it would be possible to benefit 
from the machinery developed in semiclassical physics 
(see, for example, Ref. ya and references therein) to ana- 
lyze quantum chaos in connection with multiple scatter- 
ing. As far as we know, such a semiclassical approach 
has never been considered in the context of photonic 
crystal physics but it seems to us very promising. In- 
deed, it is well-known that, due to convergence problems. 



band structure computations of metallic or semiconduct- 
ing photonic crystals are very heavy in the frequency 
range in which e{uj) < 0. The semiclassical approach 
could permit us to easily construct these band structures 
by taking into account the shortest periodic orbits involv- 
ing SP trajectories and lying into the Wigner-Seitz cell. 
Here, the properties we have found for the SP would be 
very useful. 
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